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Abstract 

We propose a simple Ginzburg-Landau free energy to describe the magnetic phase 
transition in solid 3 He. The free energy is analyzed with due consideration of the hard 
first order transitions at low magnetic fields. The resulting phase diagram contains all of 
the important features of the experimentally observed phase diagram. The free energy 
also yields a critical field at which the transition from the disordered state to the high field 
state changes from a first order to a second order one. 
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I. Introduction 

Solid 3 He undergoes a phase transition^ 1 ) at Ty = lmK into an ordered state that 
consists of ferromagnetic (100) planes and two successive planes with parallel spins followed 
by the next pair of planes where the spins are antiparallel. This state, as the magnetic 
field is increased to approximately 0.5 T, is replaced by a normal canted antiferromagnet. 
The transition to the low field uudd (up-up-down-down) phase is first order. It is also first 
order to the canted antiferromagnet phase for fields B < 0.6T. It becomes second order( 2 ) 
for B > 0.6T. The point B = 0.6T is similar to a critical point. 

The standard model^ 3 ) for solid 3 He consists of pair, triple and four spin exchanges ana- 
lyzed within a mean field theory. The model has been successful in deriving a large number 
of experimental observations. The multiple spin exchange model provides important phys- 
ical insights. It shows that because of the competing interactions (even spin permutations 
lead to antiferromagnetic and odd permutations to ferromagnetic exchange interactions, 
see for example D. J. Thouless, ref. (3)) both the Curie- Weiss constant and the bare Neel 
temperature are small. This feature for example resolves a long standing disagreement 
between exchange constants measured from T\ measurements and from Curie- Weiss tem- 
perature. While the former is an average of squares of the exchange constant, the latter 
is a linear sum. Exchange constants of different signs contribute differently to the two 
observables. It also shows that solid 3 He is a frustrated antiferromagnet. There are some 
difficulties. For example the molar volume dependence of various observables is almost 
identical, giving rise to suspicion*- 4 ) that a much smaller number of energy scales are in- 
volved. Quantum Monte Carlo calculations^ 5 ) show a weak convergence in the magnitudes 
of exchanges involving larger number of spins. Finally, multiple spin exchange models are 
sufficiently complex so that it often becomes difficult to obtain a qualitative understand- 
ing of processes in solid 3 He, something that lately has become essential as more complex 
phenomena are discovered. For example, in thermal conductivity measurements ^) there 
is evidence for magnetic defects and we hope to be able to calculate the energies of these 
defects through a Ginzburg-Landau model. 

Our aim in this paper is to propose a Ginzburg-Landau type free energy. The free 
energy is not — at least yet — derived from any microscopic Hamiltonian. Its terms are 
those allowed by the various ground states, satisfying lattice symmetries, and interactions 
between them. We study thermodynamic consequences with the aim of fine tuning its 
terms to produce a coherent, albeit phenomenological model for solid 3 He. Eventually, one 
task will be a microscopic calculation of the parameters. 
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Our starting point is a proposal by Guyer and Kumar^ 7 ) (GK). Noting that the spin 
susceptibility was close to Curie law just above the transition at low fields, they argued 
that this implied a special cancellation of the interaction energies. They further argued 
that the nature of the first order transition, hard in the sense that the entropy discontinuity 
is a major fraction of Rln2, further implied that the free energy for solid 3 He is given by 

F = Jm\ — ttiqB — TE(m) (1) 

Here B is the magnetic field, T is the temperature and E, the entropy given by 

E(£) = £n2 - h(l + S)£n(l + S) + (1 - S)£n(l - S)] (2) 

The magnetization m(r) is a vector of length fxjyS and is given by 

mix) =m + rn 2 cos(k 2 -r + n/4) (3) 

The vectors tuq -777-2 = and |m| 2 = uiq + 171?,. In our notation, when the vector sign is 
omitted, the quantity refers to the amplitude of the vector. Here k 2 is half of the reciprocal 
lattice vector in (100) direction and leads to the uudd state. The phase angle 7r/4 ensures 
that all spins are of the same length. The entropy E(S') is the entropy of a spin 1/2 system. 
In general the entropy is a consequence of proper configuration counting subject to the 
constraints on the length and the discrete nature of the spins. The expression in Eq. (2) 
reproduces the mean field equation of state. In the following we shall use the entropy as 
in Eq. (2). Note however that S is scaled so that its two possible values are S 2 = 1. The 
effect of a quadratic term in 7772 has been discussed in some detail in ref . (8) . 

In a recent preprint, ^ we have argued that a Ginzburg-Landau theory as an expansion 
in powers of the order parameter is inappropriate (and inadequate) to describe a first order 
phase transition. In a hard first order phase transition, the order parameter, at the phase 
transition jumps to a value which is easily outside the range of validity of Ginzburg-Landau 
theory. We cannot expand the entropy in Eq. (1) in powers of the order parameter and 
Eq. (1), subject to the constraint, m? = 7«q + m 2 + m 2 must be solved, if necessary, 
numerically as was already done in ref. (8). While Eq. (1) was applicable only at low fields 
and near the phase transition, the corresponding free energy for all fields is described in 
the next section. 

This paper is organized as follows: in section II, we introduce the free energy involving 
leading contributions of all of the possible ground states. Section III contains a solution 
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of the phase diagram for B = T = in the parameter space of the free energy. Section 
IV describes the results and we conclude with a summary and a discussion of possible 
improvements. 

II. The Model 

The minimal model should contain the three order parameters mo, mi, and m 2 such 
that local magnetization is given by m(r) = mo + m\cos(k\.r) + m 2 cos(k 2 .r). Here k\ and 
k 2 are respectively the wavevectors for the simple antiferromagnetic and the uudd states. 
The transition in is hard first order and we consider only terms quartic in m^- The 
order parameter m\ represents the amplitude of a simple antiferromagnetic state, present 
inthe high field part of the phase diagram. Since the phase boundary between high field 
phase and the paramagnetic phase has a positive slope > 0) in the high field phase, 
mo must have an attractive interaction with the mode m\. In general, mo, m\ and m.2 are 
all vectors. We will omit this and also any anisotropy energy terms; they are much smaller 
energy scales. We have for the free energy 

F = -Jm\ - Jim\ - J 2 m\m\ - m B - k B TZ(S) (3) 

The various = tuqSi where mo = gHjj\S\. Thus the amplitudes Sq, Si and S 2 , & t 
T = can be represented by points on a unit circle in 3 dimensions. In scaled variables 

/ = F /Jm 4 = -Si - 9l Sf - g 2 Slsl - bS - iZ(S) (4) 

Here g\ = Ji/JrriQ, g 2 = J2/J1 b = B / JrriQ and t = k^T/Jru^. The amplitude S is given 
by S 2 = S$ + Sf + Si 

At b = t = 0, we can draw a phase diagram (see Fig. 1) in the plane (gi, g 2 ) for the free 
energy in Eq. (4). At t = 0, we have S 2 = 1. For small gi and g 2 , only S 2 can be non-zero. 
Since / is monotonic in S 2 , the minimum for / must lie at S 2 = 1 so that / = — 1. We 
call this region I. In region II, g 2 is small and so we take Sq = S 2 = 0. Again since the 
minimum of / (at £1 = 1) must be /// = —g\. This minimum is lower than region I when 
g\ > 1. Thus the phase boundary I II must be g\ = 1. 

For larger values of g\ and g 2 , we expect S 2 = 0. We can write 

fm = -gi(i ~ S$) - g2S 2 (i - S^) (5) 

Its minimization with respect to Sq leads to 
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and 

,„,-^-fi(i-!) , ~<s±s£ (7) 

The boundary II III, given by fjj = fjjj is therefore g\ = g<i- The boundary I III is 
given by 

92 = 2 - gi + 2y/(l- 9l ) ~ 4 - 2(7i + 0((7 2 ) (8) 

The point {91,92) = (0,4) can also be derived from a stability analysis. The point 
(<7l> #2) = (1) 1) is a confluence of three different phases and may be called a triple point. 
These results are described in Fig. (1). 

At t = 0, the problem is still solvable, to some extent. For example, we consider the 
case of the upper critical field. In the ordered state, S2 = and S\ vanishes at the upper 
critical field. The transitions are all second order. We have (on substituting Sq — 1 — S 2 , 
and expanding in powers of S\) 



f = -b 



(91 +92)- ^ 



Sf + (g 2 + ^ Sf + O(Sf) (9) 



We see that an S\ 7^ solution exists only for b < 2(^1+^2) = b C2 - In order that the 
6 = ground state be £2 7^ 0, the ((71,(72) parameters have to come from region I. For 
(72 > 1, Eq. (8) leads to an inequality, 

6 C2 =2((7i + (7 2 )<4 v ^ (10) 

To calculate the lower critical field, b Cl , we note that the field causes the ground state 
to change from the one in region I to the one in region II. For a given set of {91,92) the 
ground state in region II contains a uniform magnetization So given by Eq. (6). In region 
II the free energy is lowered by the Zeeman interaction —6 • Sq. To order b the free energy 
in region I remains unchanged. Equating the two free energies, we get, 



b Cl = J , 292 d l- i9l t 92)2 } (11) 
! ■ ! {92-9lV 4(72 J 1 ' 

The lower critical field is negligibly small near the zero field phase boundary between 
regions I and III. It also diverges at the line g\ = (72- Experimentally the value of the 
ratio b C2 /b Cl ~ 40 and the parameters have to be rather close to the I-III phase boundary. 
While <72 < 4, g\ must be quite small. 
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To determine the critical field when the transition to the high field state becomes 
second order, we have to include the temperature t ^ 0. Since the transition is second 
order, at least for b > bt, a Ginzburg- Landau expansion of the entropy is expected to be 
valid for Si- However the magnetization Sq is large and the entropy expansion does not 
work. 

III. Phase Diagram 

We have chosen (<7i,<72) = (0-5, 1.5) for illustration. The ratio of upper critical field 
to lower critical field for this choice is 7. This is far from the experimental (see e. g. 
Godfrin and Osheroff in ref. 3) value 47. However, the qualitative features are essentially 
independent of the exact value of (gi, (72) i n a given region in Fig. (1). 

It is not possible to determine ((71,(72) from the T = critical fields alone. For a 
given ratio of B c 2/B c \, we can determine a contour in the ((71,(72) plane. If we take 
the experimental value of the B c \ = 0.45T, then b c \/t c = 0.37 using the exeprimental 
constants for the nuclear Bohr magneton. Thus b c \ = 0.53, (for t c = 1.45). Using the 
ratio for B c 2/B c \ = 47, we get b c 2 = 25. The model has an upper limit of b c 2 < 8 to 
guarantee that the ground state at T = is the uudd state. There is an easy resolution 
of this problem, namely a g shift, all 6's are replaced by a (736. This however leaves us 
unable to determine the parameters based on the phase diagram alone. In any case there 
needs to be some accomodation of the finite temperature renormalizations which in turn 
involves the mean field assumption on the entropy. We therefore are reluctant to use finite 
temperature observables. We hope to calculate temperature independent properties such 
as spin wave velocities and the energies of defects to fix the parameters at a later date. 

At b = 0, the transition occurs at t c = 1.45 into the uudd state. The entropy dis- 
continuity is almost full in accord with ref. (7). Note that if the entropy £(£) had been 
expanded in powers of S, as noted in ref. (8), the transition temperature would have been 
2.93 (t^ 1 = Yi H — 7=) and the change in S^AS^ = 1.97. In the present calculations 
AS2 = 1. Since Sq and Si are zero in the ground state, the transition temperature is 
independent of ((71,(72)- However the lower critical field does depend on these parameters. 

As the magnetic field increases, t c slowly decreases. There is a point at which all three 
phases (S2, (So, Si) and paramagnetic) coexist, a triple point. The transition between low 
field phase (S2) and high field phase (Sq, Si) is first order, as is the transition between the 
paramagnetic and high field state. These experimental features appear to be reproduced 
by the model. It came as a surprise that the model also produces a critical point. This is 
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the point (bt, tt) at which the first order phase transition between disordered paramagnetic 
phase £0 7^ 0, Si = and the high field phase (So, Si) 7^ becomes a second order phase 
transition. 

We thus see in Fig. (2) a complete phase diagram for solid 3 He magnetism. In region I, 
the ground state is the well known up-up-down-down state. In region III, the ground state 
is So 7^ 0, Si 7^ 0. Region II represents the disordered state. Again the choice of parameters 
determines the ratio b C2 /b Cl — 7. The transition at b = is independent of the parameters 
(<7l) 92) an d is at 1.4. The precise value of b Cl (at .6) clearly depends on (gi, (72)- 

In Fig. (2), the solid lines describe a first order phase transition while the dashed 
line represents a second order transition. Fig. (3) shows the discontinuity in Si(ASi) as a 
function of magnetic field along the phase boundary. We see that ASi vanishes at b = 1.25. 
For b > 1.25, the phase transition is second order. The transition temperature increases 
with magnetic field until t ~ 2.2 and b = 2.7. For b > 2.7, the transition temperature 
decreases with increasing magnetic field, reaching t = at b C2 = 2(gi + (72) = 4. 

IV. Summary 

We have shown that the phase diagram of solid 3 He can be derived from rather simple 
free energy considerations. These are extensions of a model proposed by Guyer and Kumar 
for the ordering transition at low magnetic field. The present results reproduce the essential 
features of the phase diagram including the critical point in the high field- paramagnetic 
phase boundary. 

There are two possible directions for further work, both involving the introduction of a 
space gradient dependent term in the free energy. These are (1) the derivation of magnon 
(spin wave) dispersion and (2) analysis of defects and their interaction with the magnons. 
The dynamics of spins can be easily written down using the Bloch equations for the motion 
of spins subjected to a local field derived from Eq. (4). The defects are most likely domain 
walls between metastable states frozen into the true ground state. We will return to these 
questions later. 

We acknowledge useful discussions with E. D. Adams, R. Guyer, M. Roger and N. 
S. Sullivan. This work was partially supported by the National High Magnetic Field 
Laboratory and the US Department of Energy, grant DEF G05-91-ER45462. 
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Figure Captions 

1. The phase diagram in the ((71,(72) parameter space at b = t = 0, zero magnetic field 
and temperature. In region I, the only non-zero order parameter is £2, representing 
the up-up-down-down state. All other amplitudes (Sq, Si) are zero. Similarly in region 
II only S\ 0, Sq = S2 = 0. In region III in contrast, S2 = and Sq and £1 are finite. 

2. The phase diagram in the (6, t) plane for ((71, (72) = (0-5, 1.5). The solid lines represent 
a first order phase transition while the dashed line represents a second order transition. 

3. The discontinuity in Si(AS\) at the phase transition. At small fields (6>0.5), the 
transition is clearly hard second order since AS\ ~ 0.6. However for b > 1.25, the 
transition is second order. 
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